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INTRODUCTION 


Residual  stresses,  including  macro  and  micro  stresses,  are  caused  by  the  nonuniform  plastic  flow 
of  previous  operations.  Most  manufacturing  processes  involving  thermomechanical  and  chemical 
treatment  of  components  generate  residual  stresses.  In  weapon  systems,  artificial  processes  are  used  to 
induce  advantageous  residual  stress  distribution  for  retarding  fatigue  crack  growth  and  improving  fatigue 
strength  of  the  component. 

In  recent  years,  advances  have  been  made  in  nondestructive  techniques  and  instrumentation 
employing  x-ray  diffraction,  ultrasonics,  neutron  diflraction,  and  Barkhausen  noise  techniques  for  residual 
stress  determination.  Among  these,  x-ray  diffraction  has  been  established  as  one  of  the  most  accurate  and 
rehable  methods  (refs  1-3).  The  D-IOOO-A  x-ray  stress  analyzer,  manufactured  by  Denver  X-Ray 
Instruments,  is  a  single-exposure  position-sensitive  scintillation  detection  (PSSD)  system.  It  allows  fast, 
noncontacting,  nondestructive  biaxial  surface  stress  anafysis  in  metals,  alloys,  and  other  polycrystalline 
materials.  This  report  summarizes  the  experimental  methods  employed  during  the  calibration  of  this 
instrument.  Utilizing  the  four-point  bend  apparatus  with  a  bar  specimen  made  of  ASTM  A723  gun  steel, 
x-ray  stresses  and  strain  gauge  results  were  compared,  and  the  elastic  constant  for  the  211  plane  of 
martensitic  steel  was  determined.  Single-exposure  stress  analysis  was  made  for  several  specimens.  Local 
software  development  has  been  made  to  use  the  single-exposure  PSSD  for  multiple-exposure  (sinV)  x-ray 
diffraction  stress  anatysis. 


INSTRUMENTATION 

The  D-IOOO-A  x-ray  stress  analyzer  is  a  prototype  instrument  designed  after  a  Ruud-Barrett  PSSD 
(refs  4,5).  The  instrument  utilizes  a  miniature  x-ray  tube  that  fits  inside  a  cylinder  with  a  94-mm  bore 
^ameter.  The  operating  principles  for  this  detection  qrstem,  shown  in  Figure  1,  consist  of  converting 
incident  x-rays  into  light  throu^  cadmium-zinc  sulflde  scintillation  coating;  transporting  this  light  to  a 
detector  package  through  fiber  optics  bundles;  amplifying  the  light  intensity  through  an  image  intensifier; 
and  converting  it  into  electronic  signals  through  a  reticon  photodiode  array.  The  PSSD  exploits  x-ray 
sources  of  about  100  to  200  watts  and  provides  a  shorter  data  collection  time  interval  of  any  x-ray 
diffraction  detection  ^stem. 

The  system  currently  uses  a  chromium  target  tube,  whose  characteristic  radiation  does  not 
fluoresce  any  major  alloys  in  steel.  Chromium  K-a  radiation  reflects  from  the  211  plane  of  the  body- 
centered-cubic  (BCC)  martensitic  steel  at  a  20  angle  of  156.41  degrees.  The  reticon  diode  arrays,  which 
are  solid  state  scanners,  consist  of  a  row  of  silicon  photodiodes,  each  with  an  associated  storage  capacitor 
upon  which  the  photocurrent  is  integrated.  The  device  has  sensors  on  25.4  micron  centers  with  a  sensing 
aperture  430  microns  wide.  There  are  512  diodes  in  each  of  the  two  arrays.  The  array  dimensions  are 
length  =  25.4  micron  x  512  =  1.3  cm  and  width  =  430  micron  or  0.43  mm.  The  detector,  1.3  cm  in 
length,  spans  an  angle  of  21.2  degrees  for  a  focusing  circle  of  radius  R^  =  35  mm.  To  obtain  the  required 
sensitivity,  the  silicon  diode  array  is  cooled  to  approximately  -20°C  by  a  copper  heat  sink  connected  to  a 
water  jacket.  The  cylinder  is  hermetically  sealed  and  flushed  with  nitrogen  so  that  water  vapor  does  not 
condense  on  the  cooled  parts. 

An  IBM  AT  computer  is  connected  to  the  system  for  data  acquisition,  control,  storage,  and 
anafysis.  Since  the  PSSD  produces  analog  signals,  data  have  to  go  through  an  analog-to-digital  converter 
and  associated  electronics.  The  system  software  package  uses  a  number  of  algorithms  to  correct  electronic 
and  mechanical  hardware  fluctuations,  as  well  as  x-ray  focusing  errors.  Specific  data  fitting  algorithms  are 
used  for  background  correction  and  diffraction  peak  determination. 
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SINGLE-EXPOSURE  X-RAY  DIFFRACTION  STRESS  ANALYSIS 


In  x-ray  diffraction,  stress  is  determined  indirectly  by  the  strain  it  produces.  Hooke’s  law  relates 
stress  and  strain  by  the  stiffness  and  compliance  moduli.  Stress  and  strain  are  second  rank  tensors,  while 
the  stiffness  and  compliance  moduli  are  fourth  rank  tensors.  From  symmetry  relationships  in  most  crystal 
structures,  the  number  of  independent  nonzero  constants  in  these  matrices  is  greatly  reduced.  Since  stress 
perpendicular  to  a  free  surface  is  zero,  our  interest  focuses  on  biaxial  stress  analysis. 

Characteristic  x-ray  radiations  scattered  from  the  successive  crystalline  planes  undergo  constructive 
interference  if  certain  geometric  conditions  satisfying  Bragg’s  law  are  met.  The  interplanar  lattice  d- 
spacing  acts  as  a  strain  gauge.  Diffraction  peak  location  measurements  can  be  used  to  determine  both 
residual  and  loading  stresses.  As  shown  in  Figure  2,  stress  is  determined  from  the  shift  in  diffraction 
peaks  located  on  two  opposite  sides  of  the  diffraction  cone  in  a  single  psi  tilt: 

a  =  _ M _  (1) 

(1+v)  •4/?^sm^0sin(2p) 


where  E  is  Young’s  modulus,  v  is  Poisson’s  ratio,  (Sj-Sj)  represents  the  peak  shift  between  the  right  and 
left  diffraction  peaks,  is  the  radius  of  the  focusing  circle,  /3  is  the  angle  the  incident  x-ray  makes  with 
the  specimen  normal,  and  0  is  the  Bragg  angle  of  the  unstressed  metal  (ref  6). 


MULTIPLE-EXPOSURE  X-RAY  DIFFRACTION  STRESS  ANALYSIS 

In  single-exposure  analysis  using  PSSD,  a  single  psi  tilt  gives  12  points  on  the  sinV  plot.  Multiple- 
exposure  x-ray  diffraction  stress  analysis  depends  on  the  measurements  from  multiple  psi  angles  of 
inclination.  Using  the  PSSD,  N  number  of  psi  tilts  gives  2N  number  of  points  on  the  sin^  plot  for  stress 
determination  by  the  least  squares  fit  of  a  line.  Multiple-exposure  stress  analysis  using  PSSD  has 
advantages  over  both  the  single-exposure  technique  using  PSSD  and  the  conventional  multiple-exposure 
technique  using  a  single  detector.  Local  software  development  has  been  made  to  use  the  PSSD  for 
multiple-exposure  biaxial  stress  analysis. 


EXPERIMENTAL  METHODS 
Surface  Preparation 

X-ray  penetration  depends  on  the  x-ray  wavelength  and  the  absorption  coefficient  of  the  specimen. 
For  chromium  x-rays  on  steel  samples,  penetration  is  about  twenty  microns.  To  make  an  accurate  residual 
stress  determination,  removal  of  the  specimen  surface  affected  by  cold  working,  e.g.,  grinding  and 
polishing,  corrosion,  oxidation,  and  grease  and  polymer  coatings,  is  necessary.  This  is  done  by  machine 
and  hand  polishing,  metallographic  pohshing,  and  then  electropohshing  of  the  specimen  surface.  We 
experimented  with  current  density  settings,  concentrations  of  the  polishing  bath,  and  used  both  lead  and 
stainless  steel  mesh  as  the  cathode  in  the  electropolishing  procedure.  Figure  3  shows  the  electropolishing 
setup  used  in  our  experiments. 
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Detector  Gain 


When  a  constant  x-ray  flux  is  incident  on  the  detector,  pixel-to-pixel  differences  in  total  detection 
efficiency  give  rise  to  gain  variations.  The  variations  in  the  scintillation  coating  thickness,  scintillation 
conversion  efficiency,  variable  amplification  of  the  image  intensifier,  nonuniform  sensitivities  of  the  silicon 
diode,  fiber  optics  variations,  and  coupling  factors  can  contribute  to  the  nonuniformity  of  the  gain  signals. 

To  make  pixel-to-pixel  efficiency  correction,  gain  signals  were  obtained  by  irradiating  a  glass 
surface  with  the  x-ray  beam  at  the  same  geometiy  the  stress  measurements  were  to  be  made.  Since  glass 
is  an  amorphous  solid,  it  should  have  a  uniform  diffraction  intensity  distribution  versus  angle  26  in  the  far 
back  reflection  range. 

To  investigate  the  uniformity  and  stability  of  the  gain  signal,  a  log  of  gain  signal  intensity  versus 
diode  array  element  was  plotted  at  various  integration  times  as  shown  in  Figure  4  and  at  various  tube 
currents  as  shown  in  Figure  5.  Large  pixel-to-pixel  intensity  variations  were  observed,  especially  in  the 
right  detector.  As  shown  in  the  following  sections  of  this  report,  diffraction  profiles  have  to  be  corrected 
for  nonuniform  gain  variations.  Figure  6  shows  that  electronic  drift  in  the  detectors  causes  a  gain  shift. 
Two  traces  of  gain  signals  versus  diode  array  elements  were  recorded  at  subsequent  times  using  the  same 
integration  time,  tube  current,  and  voltage,  and  they  were  plotted  on  the  same  graph.  Figure  7  shows  the 
effect  of  the  noisy  gain  on  the  noise-  and  gain-corrected  diffraction  profiles.  Gain  stability  and  statistics 
are  essential  in  accurate  peak  location  determination.  In  this  work,  gain  signals  were  obtained  in  5-  to  10- 
second  integration  times.  It  was  continuously  and  carefully  monitored  for  possible  electronic  gain  drift. 

Detector  Noise 


Possible  hght  leaks  and  the  electronic  dark  current  both  in  the  image  intensifier  and  diode  arrays 
contribute  to  the  detector  noise,  which  is  detected  when  no  x-rays  are  applied.  The  dark  current  in  the 
diode’s  arrays  increases  with  the  temperature  of  the  arrays.  Detector  noise  signals  were  collected  for  both 
the  left  and  right  detectors  in  the  PSSD  and  were  used  to  correct  for  the  gain  and  profile  readings. 

Figure  8  is  a  plot  of  the  average  noise  variation  with  detector  integration  time  for  the  left  and 
right  detectors  in  the  PSSD.  This  plot  indicates  that  the  right  detector  has  a  noise  offset  of  87,  and  the 
left  detector  has  a  noise  offset  of  25.  To  equalize  the  noise  offsets,  potentiometers  within  the  instrument 
have  to  be  adjusted.  At  an  integration  time  near  forty  seconds,  the  curve  begins  to  show  an  exponential 
behavior.  At  an  integration  time  less  than  ten  seconds,  which  was  used  throughout  this  work,  a  linear 
relationship  between  noise  and  the  integration  time  is  observed.  A  formula  was  derived  from  this  plot  and 
system  software  was  modified  to  normalize  integration  time  differences  for  the  noise,  gain,  and  profile 

Noise(t^  =  Noise(t^  +  •  Slope  (2) 


where  tj  and  tj  are  two  different  integration  times. 

Diffraction  Profile 


The  intensity  of  the  diffraction  profile  is  proportional  to  the  intensity  of  the  characteristic 
radiation,  which  is  in  turn  linearly  dependent  on  the  tube  current  and  nonlinearly  on  the  tube  high  voltage. 
In  this  work,  a  tube  voltage  of  80  kilovolts  and  a  current  of  6  milliamps  were  used.  Integration  time  for 
the  diffraction  profiles  ranged  between  2  and  10  seconds.  Enough  sampling  time  should  be  allowed 
without  saturating  the  charge  capacitor. 
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Diffraction  peak  profiles  are  noise-corrected  and  divided  by  the  gain  efficiency  to  normalize  the 
pixel-to-pixel  gain  sensitivities.  System  software  was  modified  to  perform  this  operation. 

P(i)  =  (Profile(i)  -  Noise(i))KGam(i)  -  Noise(i))  (3) 

To  obtain  accurate  stress  measurement,  sufficient  integration  time  should  be  allowed  for  the  profile  and 
the  gain  signal.  Diffraction  peak  profile  detenmnation  was  found  to  be  difficult  for  noisy  gain. 

Specimen  Displacement  Error 

If  the  irradiated  diffraction  volume  is  displaced  from  the  center  of  the  diffractometer,  there  is  a 
relative  peak  shift  between  psi  =  0  and  psi  =  ^  due  to  specimen  displacement  error.  The  Denver  x-ray 
stress  anal3^er  has  a  fl-goniometer  geometry,  and  the  error  in  26  location  is  given  by  Eq.  (4): 

6(A2e)  =  (A20)^^  -  (A20)^^^  =  f^VA;O-cos0 


sin0 


/&m(0+tlf) 


(4) 


where 

R  =  R  cos(t|t+(90°  -0)) 
"  cos(ilf-(9O°-0)) 


In  Eq.  (4),  is  the  diffractometer  radius,  AX  is  the  specimen  displacement  error,  6  is  Bragg’s  angle,  and 
ijj  is  the  psi  tilt.  When  the  specimen  is  too  high,  AX  is  positive.  When  the  specimen  is  too  low,  AX  is 
negative.  In  the  sin^^  plot,  high  specimen  position  contributes  to  negative  slope  and  compressive  stress 
state.  Low  specimen  position  contributes  to  positive  slope  or  tensile  stress  state.  Precision  specimen 
alignment  by  a  micrometer  dial-gauge  indicator  or  laser-sensor  indicator  is  suggested.  ASTM  recommends 
zero  stress  standard  of  loose  powered  iron.  In  this  work,  stress-free  400  mesh  iron  powder  was  used  to 
obtain  zero  stress  correction  and  to  check  for  specimen  alignment.  Diffraction  peak  profiles  for  zero 
stress  powder  are  given  in  Figures  9  and  10. 

Peak  Location  Determination 

Before  the  diffraction  peaks  can  be  located,  the  raw  diffraction  peak  intensity  profiles  must  be 
processed  by  several  correction  algorithms  as  provided  by  the  PSSD  system  software  (ref  7): 

1.  Incoherent  bremsstrahlung  background  radiation  subtraction. 

2.  Lorentz  polarization  absorption  intensity  corrections. 

3.  Planer  instead  of  curved  fiber  optics  surfaces  were  used.  To  correct  for  the  fiber  optics 
alignment  in  the  location  of  the  center  of  the  x-ray  sensitive  fiber  optics,  a  semi-empirical 
parabolic  regression  method  was  used. 

4.  Iterative  correction  for  the  uncertainty  in  the  detector  circle  to  specimen  distance  R^,. 
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Li  and  Sun  (ref  8)  studied  contributions  to  the  background  radiation  in  x-ray  stress  analysis 
following  derivations  from  diffraction  physics.  They  proposed  a  nonlinear  equation  for  background 
correction  dependent  on  20  and  the  psi  tUt,  which  agreed  well  with  the  experimental  background  curve. 
Furthermore,  we  propose  that  the  cross  correlation  technique  be  used  to  determine  the  diffraction  peak 
shift  for  weak  and  noisy  diffraction  profiles  that  might  occur  in  certain  crystalline  structures.  In  this  work, 
a  straight  line  background  fitting  algorithm  was  used.  The  user  determines  a  range  for  the  starting  and 
ending  pixel  locations  to  make  the  background  correction.  The  user  can  select  the  percentage  of  the 
upper  peak  intensity  profile  in  the  analysis  for  the  parabolic  fit.  We  found  that  70  to  80  percent  of  the 
maximum  peak  intensity  is  good  for  the  parabolic  peak  determination. 

Cross  Correlation  Technique 

The  cross  correlation  technique  can  be  applied  to  determine  diffraction  peak  shift,  especialfy  for 
weak  and  noisy  diffiraction  profiles  that  might  occur  in  certain  crystalline  structures.  By  cross  correlating 
diffraction  peak  intensity  versus  diode  array  element  signal  at  psi  =  0  with  signals  at  psi  =  and  psi  = 
tp2,  respective^,  residual  stress  can  be  calculated  fi-om  the  diffraction  peak  shift.  An  alternate  method  is  to 
cross  correlate  the  left  detector  signal  (Fj)  and  the  right  detector  signal  (Fj).  The  cross  correlation 
function  for  each  ip  tilt  is  given  by 


=  E  A 


512 


i=l 


(5) 


where  i,j  are  diode  array  element  indices.  From  the  location  of  the  maximum  in  the  cross  correlation 
function,  peak  shift  (Sj-Sj)  in  Eq.  (1)  can  be  obtained. 


FOUR.POINT  BEND  CALIBRATION  AND  ELASTIC  CONSTANT  DETERMINATION 

A  four-point  bend  instrument  was  used  to  compare  the  x-ray  and  strain  gauge  stress 
measurements.  A  strain  gauge  was  attached  to  a  6  by  1  inch  ferritic  steel  specimen  made  fi-om  ASTM 
A723  gun  steel,  and  both  tensile  and  compressive  stresses  were  applied.  In  Figures  11  and  12,  the  raw  and 
corrected  peak  profiles  for  the  ferritic  steel  specimen  are  shown  for  the  right  and  left  detectors, 
respectively.  In  Figure  13,  the  x-ray  stress  is  plotted  versus  strain  gauge  stress.  The  elastic  stress  constant 
was  determined  from  the  slope  of  the  curve  as  follows: 

-  _  0.0254 

V  1+v  j  ^4R^sin^0sm(2p)^ 


/  N 

a 

1+v 

y,  2  i) 

where  K  = 


0.0254 

4R^sin^0sin(2P) 


(6) 


Our  experiment  indicated  that  the  constant  Eijj,/(l-»-w)  =  23.9  •  10®  psi.  The  quantity  was  determined  to 
be  23.5  •  10®  psi  for  the  same  specimen  using  the  stress  analyzer  at  the  Materials  Research  Laboratory, 
Penn^lvania  State  University  during  our  visit. 
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STRESS  ANALYSIS  OF  M106  SHELL  SAMPLE 


In  Figure  14,  the  raw  and  corrected  right  detector  peak  profiles  for  the  inside  surface  axial  stress 
of  an  M106  shell  section  are  shown.  The  left  detector  inside  surface  hoop  stress  profiles  of  an  M106  shell 
section  are  shown  in  Figure  15.  In  Figure  16,  the  multiple-exposure  technique  was  used  for  the  M106 
shell  inside  axial  stress  determination.  Psi  tilts  at  12,  20,  25,  and  30  degrees  give  pairs  of  data  points  (1,5), 
(2,6),  (3,7),  and  (4,8)  on  the  graph.  D  is  the  interplanar  spacing  of  the  planes  at  nonzero  psi  and  D,,  is  the 
interplanar  spacing  of  an  unstressed  metal.  The  value  was  obtained  from  the  literature  (ref  9).  Linear 
regression  analysis  was  used  to  obtain  the  slope  for  residual  stress  determination.  In  Table  1,  the  stress 
results  measured  at  Benet  Laboratories  are  given  along  with  measurements  from  Pennsylvania  State 
University.  Errors  quoted  represent  data  dispersion  in  successive  measurements,  including  statistical  error, 
systematic  error,  and  error  due  to  goodness  of  fit  in  the  multiple-exposure  analysis.  Residual  stress 
accuracies  are  expected  to  be  approximately  ±  10  Ksi.  The  measurements  are  in  good  agreement— Benet’s 
values  were  slightly  higher  than  Pennsylvania  State’s. 

Table  1.  Residual  Stresses  in  a  Shell  Specimen 


Stress 

Orientation 

Bea6t  Labs* 

Penn  State** 

Bea&t  Labs*** 

Stress  (Ksi) 

&ress  (Ksi) 

Stress  (Ksi) 

inside  hoop 

-42.1  ±  3.3 

-39.1  ±  3.1 

43.5  ±  2.6 

inside  axial 

43.4  ±  3.3 

34.6  ±  1.8 

-40.3  ±  3.9 

outside  hoop 

-39.9  ±  1.4 

-35.0  ±  1.6 

46.3  ±  3.2 

outside  axial 

-39.1  ±  3.2 

-35.3  ±  1.7 

43.4  ±  2.6 

*Denver  D-IOOO-A  Single-Exposure  Stress  Analyzer  (local  software  modification) 
**Pennsylvania  State  University  Single -Exposure  Stress  Analyzer 
***Denver  D-IOOO-A  Multiple-Exposure  Measurements  (local  software  development) 


CONCLUSION 

Experimental  methods  using  Denver  D-IOOO-A  single-exposure  position-sensitive  scintillation 
detection  system  (PSSD)  are  given  in  this  report.  Characteristics  of  noise,  gain,  and  diffraction  profiles 
were  studied  as  a  function  of  system  parameters  such  as  integration  time,  number  of  iterations,  and  x-ray 
tube  current  and  voltage.  A  four-point  bend  experiment  was  performed  for  an  ASTM  A723  steel 
specimen.  Elastic  constant  of  the  211  plane  of  BCC  steel  was  determined.  Local  software  development 
was  made  to  allow  the  single-exposure  PSSD  to  run  in  a  multiple-exposure  mode.  The  technique  improves 
accuracy  in  the  sinV  biaxial  stress  analysis  using  the  PSSD.  Residual  stress  analyses  were  made  for  several 
ordnance  steel  specimens. 

Suggested  software  enhancements  include  cross  correlation  technique  for  peak  shift  determination; 
capability  to  determine  peak  position  by  profile  fit  to  Gaussian,  Lorentzian,  and  Pearson  VII  profile  types; 
layer  removal  stress  corrections;  and  triaxial  stress  analysis  for  complete  stress  tensor  determination. 
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Suggested  system  hardware  improvements  include:  (1)  uniform  and  stable  gain  by  recoating  the 
scintillation  detector;  (2)  stable  current  and  voltage  power  supply  and  associated  electronics; 

(3)  precision  specimen  alignment  by  micrometer  dial-gauge  or  laser-sensor  indicators;  (4)  automated  psi 
motion  control  for  multiple-exposure  stress  analysis;  and  (5)  automated  radiation  shutter.  By  the 
application  of  an  extended  or  mobile  photodiode  array  and  x-ray  tubes  of  other  wavelengths,  Denver  D- 
1000-A  PSSD  can  be  expanded  to  study  other  metals,  alloys,  and  composite  materials. 
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Figure  1.  Schematic  of  the  D-IOOO-A  position  sensitive  scintillation  detection  system. 
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Figure  3.  Drawing  of  electropolishing  apparatus. 
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Figure  4.  Left  &  right  detector  gain  for  various  integration  times. 
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Figure  5.  Left  &  right  detector  gain  for  various  tube  currents. 


20  (Diode  Array  Element) 


Figure  6.  Left  &  right  detector  gain  signals  indicating  shift  in  peak  locations  due  to  electronic  drift. 
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Figure  7.  Effect  of  noisy  gain 
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Figure  9.  Right  detector  diffraction  proffles  for  zero  stress  400  mesh  iron  powder. 
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Figure  10.  Left  detector  diffraction  profiles  for  zero  stress  400  mesh  iron  powder. 
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Figure  11.  Right  detector  diffraction  profiles  for  ferritic  four-point  bend  specimen. 
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Figure  12.  Left  detector  diffraction  profiles  for  ferritic  four-point  bend  specimen. 
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Figure  13.  Four-point  bend  verification  of  x-ray  stress  versus  strain  gauge  stress. 
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Figure  14.  M106  shell  section  inside  surface  axial  stress  right  detector  diffraction  profiles. 
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Figure  15.  M106  shell  section  inside  surface  hoop  stress  left  detector  diffraction  profiles. 
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Figure  16.  M106  shell  section 
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